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Abstract 

We investigate the possible measurement of parity-violating spin asymmetries in jet pair pro- 
duction in proton-proton collisions at the Brookhaven Relativistic Heavy Ion Collider (RHIC), 
with the goal to constrain longitudinally polarized quark and antiquark distribution functions. A 
measurable asymmetry could be observed in the vicinity of the massive weak W boson resonance, 
where the parity-violating signal appears above the parity-conserving background and is enhanced 
by interference of the strong and electroweak production amplitudes. We discuss the potential 
for such measurements, perhaps the first opportunity to measure a parity-violating asymmetry in 
proton-proton collisions at RHIC. Sensitivity of this measurement to the polarization of down-type 
antiquarks is demonstrated. 

PACS numbers: 12.38.Bx, 13.87.-a, 13.88.+e 
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I. INTRODUCTION 



The asymmetry in the rapidity-dependent cross section for weak (W) boson production in 
unpolarized high-energy pp collisions is a valuable instrument [H for measuring the Bjorken x 
dependence of the ratio u(x, M w )/d(x, M w ) of the up-quark and down-quark quark parton 
distribution functions (PDFs) in a proton at the very hard scale set by the weak boson mass 
Mw- Demonstration of the utility of this prediction awaited sufficient luminosity [1, 0], 
but it was not long before the measured rapidity asymmetries were used routinely in global 
analyses as important constraints on unpolarized PDFs 0, 0|. The method relies on the 
accepted Drell-Yan mechanism by which the W ± boson is produced in the standard model, 
e.g., principally u + d — > W + and d + u — > W~ at leading order in quantum chromodynamics 
(QCD) perturbation theory. Thus, at least qualitatively, a W + moving with a large rapidity 
y is produced from projectile u quarks carrying large x, whereas a W~ moving with a 
large y is produced from projectile d quarks carrying large x. The ratio of the unpolarized 
W + and W~ cross sections at large y thus provides important constraints on the ratio 
u(x, M w )/d(x, M w ) at x between 0.1 and 1. 

Extension of similar arguments to polarized pp scattering 

a, a a a a QjJ adds a 

promise of direct access to the longitudinal spin-dependent quark PDFs Aq(x, M w ) via the 
(V — A) coupling of the quark constituents to the W. Only one of the incident protons need 
be longitudinally polarized in order for this parity-violating measurement to be effective. 
In this paper we address W production in longitudinally polarized proton collisions at the 



Brookhaven Relativistic Heavy Ion Collider (RHIC) at collision energy ^fs = 500 GeV [12 
We focus on the dijet decay mode of the W, pp — > (W — > jj) + X, and compute fully 
differential cross sections in jet pair invariant mass Q, jet pair rapidity y, and jet transverse 
momentum pxj for all signal and background processes. 

The parity- violating dijet production cross section, Act = cr(p^p) — cr(p^p), where 
represents a proton with its spin aligned with (against) the proton's direction of motion, 
is dominated by the ud+ du — > W + process. The single-spin asymmetry 

A = a(p^p) - (?{p^p) ,^ 
cr(p^p) + (j(p^p) ' 

(as a function of the rapidity y ) is sensitive to Au(x, Mw) at x — > 1 at large positive y 
values, and to Ad(x, M w ) at x — » M w /s at large negative y values, as we demonstrate 
quantitatively. Both Am and Ad PDFs are not known well in the quoted x regions, as they 
are not constrained by the existing data (mostly from polarized deep inelastic scattering). 
The dijet asymmetry Al may probe these PDFs with a relatively small event sample. 

Commonly one detects a W boson by observing a leptonic decay W — > £ + u, with its 
characteristic Jacobian peak in the transverse momentum pr of the charged leptons £ = e 
or fi, and the missing transverse energy carried by the invisible neutrino v. 1 The lepton 
decay mode is observable at all colliders, including RHIC 1 1 61 ] . Such signals are clean 



and impose tight constraints on the PDFs, but the small decay branching fraction penalty 
is felt keenly in the event rate. The luminosity in polarized proton scattering tends to be 
less than in the unpolarized case, and the W production cross section at the envisioned 
RHIC center-of-mass energy is not as great as at a higher-energy collider. Correspondingly, 



1 For a general discussion of production and observation of W bosons at hadron colliders, see, e.g., Ref. 
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lepton decay events come at a premium and most likely would require several years of RHIC 
running. With this in mind, it is attractive to investigate the a priori disfavored alternative 
that one could observe W production in the hadronic decay mode W — > 2 jets precisely 
the topic of this paper. 

An advantage of the hadronic W decay compared to the leptonic decay is a larger cross 
section, enhanced by Br(W — > hadrons)/Br(W — > Iv) « 6. Another advantage is the possi- 
bility to directly determine the invariant mass and rapidity of the W boson, approximately 
equal in this case to the invariant mass (Q) and rapidity (y) of the high-p^ jet pair. The 
boson-level y distribution does not suffer from the smearing of the PDF dependence in the 
charged-lepton distributions |l4|] caused by spin correlations in the course of the W boson 
decay. 

The hadronic decay mode of the W boson has an admitted drawback in that background 
contributions to jet pair production are substantial. Hadronic jets are produced most abun- 
dantly through QCD and electromagnetic (virtual photon) interactions, rather than through 
W interactions, especially via the scattering involving real gluons, which grows rapidly with 
collider energy. The backgrounds present a major obstacle for observing unpolarized elec- 
troweak (W and Z mediated) jet production at the Fermilab Tevatron pp collider and the 
CERN Large Hadron Collider [3, [l^|. At the intermediate energy of RHIC, the back- 
ground/signal ratio (25-30 or more) is comparable to that at the CERN pp Super Proton 
Synchrotron (SPS) with the energy y^i = 630 GeV, where a 3a enhancement from W and 
Z bosons in the dijet mass (Q) distribution was observed 20| with very low integrated 
luminosity (C = 0.73 pb _1 ). 

We argue that the main backgrounds can potentially be controlled at RHIC because of 
their smooth (non-resonant) kinematical dependence and symmetry with respect to spatial 
reflections. Due to the smooth behavior, side-band subtraction can be used to advantage 
to estimate the unpolarized background accurately in the W signal region by extrapolation 
from the side bands of the Q distribution. The parity symmetry of the dominant background 
leads to its cancellation in the single-spin cross section, up to a small false asymmetry due 
to the uncertainty in the proton beam polarization. Thus the parity- violating Al preserves 
the clear pattern of PDF flavor dependence typical of W boson production and allows us 
access to Au(x, Mw) and to Ad(x, Mw) in well-defined kinematic regions. 

It is not possible to distinguish W + and W~ contributions in the jet pair decay mode. 
This apparent limitation does not disfavor the hadron decay mode with respect to the 
leptonic decay mode in situations in which charged lepton identification is not available, 
such as at the earlier stages of RHIC spin program. Full separation of the W and Z signals 
is not possible as a consequence of limited jet invariant mass resolution (of order 5-10 GeV). 
However, the Z contribution, with its relatively small event rate, has only mild effect on the 
PDF dependence of jet pair cross sections. 

We evaluate all scattering processes for production of jet pairs (approximated by a pair of 
final-state partons) at the leading order in the QCD and electroweak coupling strengths, a s 
and o>ew- Spin-dependent distributions in dijet mass (Q) and jet transverse momentum (prj) 
at this order are reported in Ref. [171]. In our paper, we focus on fully differential distributions 
and comparison of dijet production channels, issues essential for the observation of At. We 



compute these distributions using a modified Monte-Carlo integrator MadEvent [2jJ, |22|, [23 



The presence of both strong and electroweak contributions to the same final state means 
that their constructive and destructive interference is possible. We compute the magnitude 
of the interference effect in both the unpolarized and the single-spin cross sections. Although 
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the interference effect enhances the single-spin cross section considerably, it does not alter 
the sensitivity of Al to the polarized PDFs. 



In a recent paper [24|, Arnold et al. also consider jet production in polarized proton 



proton collisions at RHIC energies. Our focus and emphasis are somewhat different from 
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theirs. In our work, we compute fully differential cross sections of two jets, whereas Ref. 
investigates the single-jet inclusive double-differential cross section da / '(dpxjdyj), integrated 
over the full range of dijet mass Q. All background contributions are included in the un- 



polarized denominator of A L in Ref. 24], reducing predicted asymmetry values to about 
1% or less. In our study, we find that selection cuts on Q and other kinematical variables, 
combined with the subtraction of the background in the signal region 70 < Q < 90 GeV, 
are very effective in suppressing unpolarized dijet production. Owing to the use of these 
techniques our predictions for spin asymmetries are larger at least by an order of magnitude, 
and our appraisal of the jet pair mode is more optimistic. 

A tag on final-state charmed hadrons is proposed in Ref. H| as a means to improve the 
signal over the background. In principle charm tagging eliminates the dominant gluon 
and light-quark parity-conserving processes, while preserving a large fraction [of order 
Br (W — > cX) /Br (W — > all) ~ 33%] of electroweak signal events. In practice, the util- 
ity of this approach is diminished by the low efficiency of the experimental identification 
of final-state charmed hadrons. Charm scattering contributions are included in our Monte- 
Carlo calculation, so that charm tagging can be explored together with the other handles at 
the fully differential level. 

In Sec. HIl we discuss the basic QCD processes that produce a pair of jets. We show figures 
that illustrate the size of various QCD and electroweak contributions to the cross section, 
as a function of the dijet invariant mass and other variables. To obtain our quantitative 
estimates of polarization asymmetries in Sec. IHH we must adopt parametrizations of spin 
dependent parton distributions. We choose to work with the de Florian-Navarro-Sassot 
parametrization [25|. Using these PDFs, we explore sensitivity of the predictions for A L to 
variations in the spin-dependent PDF parametrizations. 

Conclusions are summarized in Sec. IIVI Our quantitative predictions show that the 
predicted magnitude of Al is sufficiently large relative to the anticipated uncertainties, so 
that significant measurements should be possible, with good prospects for discrimination 
among different polarized PDFs. We urge experimental study of Ai(y) in hadronic decays 
of W at RHIC. 



II. UNPOLARIZED JET PAIR PRODUCTION 
A. General remarks 

We consider inclusive jet pair production in proton-proton collisions at = 500 GeV, 
VP ~^ jjX, approximated by 2 — > 2 exchanges of the standard model bosons (V = g, 7*, 
W ± , and Z°) in the s, t, and u channels, including the interference between the amplitudes 
with different types of the bosons. The unpolarized and single-spin cross sections, a and 
Act, for production of an energetic jet pair are given by 

da _ *!££*L MxuQ)Mx , iQ)i (2) 



dftidp 2 a fi d dpidp 2 
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and 

where pi and p 2 are the momenta of the jets, Q 2 = (pi + P2) 2 , £1,2 = Qe ±v l\fs-, V = 
In l(Q° + Q 3 ) I (Q° — Q 3 )] /2, (A)a(ab — >• cd) is the perturbative cross section in the ab — > erf 
parton scattering channel, and f a /p(x,Q) (Af a / p (x,Q)) is the unpolarized (longitudinally 
polarized) PDF. The QCD renormalization and factorization scales are set equal to the 
invariant mass Q of the dijet. The z axis follows the direction of the polarized proton beam. 
In some cases, we use the notation cr(V) for the partial cross section that involves the same 
boson species V in the scattering amplitude and its conjugate, without the interference with 
other types of bosons. 

The lowest-order cross sections for s-channel W + (W~) production via qiq^ — > W ± — > 
(the dominant parity-violating contribution) are 

da(W, s-channel) ^ r . _ . ^. . . lS 

- = l^Cij Qi{xi,Q)q j {x2,Q)+q j {xi,Q)qi{x 2 ,Q) (4) 



dQ 2 dy 



1 -j 



and 



L nel) = ^ c ^ \ Aq ,^ ( X2j Q) + Ag^xx, Q)^(z 2 , <2)1 , (5) 



dQ 2 dy 



where 



7r Q 



= l± ^W * ,y X2 y ,y 



125 (Qa-M^^ + r^/M, 



VK k,l 



otEW = Qwl^ 11 i s the electroweak coupling strength, and IV are the boson mass and 
width, and Vij = Vji is the quark mixing (Cabibbo-Kobayshi-Maskawa) matrix. The indices 
i, j, k, I are summed over flavors of quark-antiquark pairs carrying the net electric charge 
of the W + or W~ boson. At the energy of RHIC, where contributions with initial-state s, 
c, b (anti-)quarks are much smaller than those with u and d (anti-)quarks, the W + and W~ 
cross sections are well approximated by 



da(W + , s-channel) 

= ^ud 



dQ 2 dy 

da(W~, s-channel) 
dQ 2 dy 
dAa(W + , s-channel) 
dQ 2 dy 

and 

dAa(W~, s-channel) 
dQ 2 dy 



c 



u(xi, Q)d(x 2 , Q) + d(x\, Q)u(x2, Q) + small terms; (6) 

C ud [d(xi, Q)u(x 2 , Q) + u(xi, Q)d(x 2 , Q)\ + small terms; (7) 
-Au(xx, Q)d(x 2l Q) + Ad{x x , Q)u(x 2 , Q) 



ud 



small terms; 

(8) 



C u d[—Ad(x 1 ,Q)u(x2,Q) + Au(xx,Q)d(x 2 ,Q)] + small terms. (9) 



These resonant cross sections are combined with non-resonant (t— and w-channel) contribu- 
tions to form a(W), the full W-boson contribution of order o? EW . In a similar manner, we 
construct pure Z and 7* contributions, also of order a 2 EW . In addition, one must include 
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the cross section for the interference between the electroweak bosons (of order a EW ), and 
between the electroweak bosons and gluons (of order otsotEw)- Finally, there is a pure QCD 
cross section, of order a 2 . Explicit matrix elements for all these processes can be found in 



Ref. 17 



Our fully differential spin-dependent cross sections for 2 — > 2 jet production are com- 



21, 22, 23 



puted in all channels with the help of the programs MadGraph and MadEvent 
MadGraph is a program for automatic generation of tree-level cross sections in the standard 
model and its common extensions. MadEvent realizes phase-space integration of these cross 
sections. Internally these programs operate with helicity-dependent scattering amplitudes 



obtained using the HELAS library [26|. In a typical setting, the amplitudes are summed 
over all helicity combinations to produce spin-averaged cross sections. We modified MadE- 
vent to also allow evaluation of single-spin cross sections for an arbitrary scattering process, 
including jet pair production. 2 This modified program is employed to evaluate both the 
numerator and denominator of Al- 

In hard-scattering cross sections, we include contributions from four quark flavors (u, d, s, 
and c), in accordance with the default MadGraph choice. The unpolarized and longitudinally 
polarized PDFs are taken from the Martin-Roberts-Stirling-Thorne (MRST'2002 NLO Q) 



and de Florian-Navarro-Sassot (DNS'2005 [25||) sets, respectively. The MRST'2002 NLO 



PDFs are chosen because they satisfy positivity conditions with the DNS'2005 PDFs. 



B. Kinematic distributions of unpolarized cross sections 

Spin-averaged jet production is dominated by the continuous event distribution from 
QCD and electromagnetic scattering (involving only g and 7*). These involve, in order of 
their magnitudes at 60 < Q < 100 GeV, qg — * qg, gg — > gg, qq' — > qq', qq — > qq, and 
smaller scattering contributions, where q stands for both quarks and antiquarks, and qq 
(qq') stands for scattering of the same (different) quark flavors. The two largest jet cross 
sections produced by qg and gg scattering cancel when a cross section difference is taken to 
compute the single-spin asymmetry. Parity violation needed to obtain a non-zero numerator 
of Al in Eq. (pp) arises solely from qq contributions with intermediate W and Z bosons. The 
s-channel parity-violating scattering amplitudes are enhanced resonantly when the dijet 
invariant mass Q is close to M w (Mz), in the Q range we will call "the signal region". 
Even in this region, the spin-averaged W and Z contribution constitutes at most a few 
percent of the full event rate, suggesting that the most straightforward measurement (not 
separating the signal and background contributions) would result in small spin asymmetries. 
The magnitude of Al may be enhanced by applying a "side-band subtraction" technique, 
i.e., by measuring the large parity-conserving background at values of Q outside of the signal 
region and subtracting it from the denominator of Al inside the signal region. 

Let us now turn to the figures illustrating these observations. The spin-averaged differ- 
ential cross sections for various combinations of scattering channels are shown in Fig. [TJ We 
focus on the distributions in the dijet invariant mass Q, the transverse momentum prj of 
the jet, and the difference Ay of jet rapidities, Ay = jyi — 2/2] 5 shown in Figs.[T](a,b,c). Cross 
sections denoted "W(a EW ) v , u Z(a Ew y\ and "7*(a|;vi/)" represent the pure W boson contri- 



2 The modified MadEvent program for the calculation of polarized cross sections is available from one of 
the authors (P.M.N.) by request. 
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Figure 1: Unpolarized cross sections for pp — > jjX at y^s = 500 GeV in various scattering channels, 
plotted as a function of (a) dijet invariant mass Q, (b) jet transverse momentum prj, and (c) 
difference of jet rapidities Ay = \yi — j/2 1 , for the cuts specified in the figures. 
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bution, the Z boson contribution, and the 7* contribution, without the interference terms. 
The sum of all electroweak and interference cross sections is shown by a lower solid line la- 
beled "a s acEw + a Ew"- The full cross section, which also includes the pure QCD contribution 
of order a 2 , is shown by the upper solid line with the label "a| + a s aEw + a Ew"- 

In most figures, we impose constraints prj > 25 GeV and \yj\ < 2 to reproduce ap- 
proximately the acceptance of the STAR detector [1^|. 3 As seen in Figs. [Qa) and (b), 
W (Z) production receives a resonant enhancement at Q ~ Mw (Mz) and prj ~ Mw/2 
(M z /2). The QCD and electromagnetic backgrounds fall smoothly with both Q and pxj- 
To focus on the parity- violating asymmetry, we impose a selection 70 < Q < Qmax-, with 
Qmax = 90 — 100 GeV, in Figs.QJb) and (c). A value of Q m ax below Mz is generally preferred 
in order to suppress the Z cross section (constituting about 1/6 of a(W) at their respective 
mass poles according to Fig. [H(a)) and emphasize the PDF dependence typical for the W 
boson contribution. 

In 2 — > 2 scattering, Ay is related to the scattering angle 9* in the jet pair rest frame as 
tanh(Ay/2) = cos#*. Signal and background processes are characterized by different spin 
correlations between the initial- and final-state particles and, therefore, different shapes of 
da I <i(cos #*) [or da/dAy]. The kinematic differences between the signal and background 
da/dAy distributions are traced largely to real-gluon emissions present in the QCD back- 
ground process, but not in the electroweak processes (see below). Figure [T](c) shows that the 
full 0(a 2 s + a s aEw + ®ew) cross section is peaked strongly at large \Ay\ , while the signal 
0{a 2 EW ) and 0(a EW + a s a E w) cross sections have flatter \ Ay\ dependence. In the dijet rest 
frame, the 0(a 2 s + a s aEw + a Ew) cross section peaks more strongly at cos6>* — > ±1 as 
compared to the 0(a 2 EW + a s a>Ew) cross section. 

The full electroweak+interference cross section constitutes 2-3% of the full cross section 
at I Ay I = 0, but it drops below 1% at \Ay\ « 2. This difference could be exploited to 
enhance the signal/background ratio. For simplicity, we don't impose this selection in the 
rest of the paper. 

In Fig. [2] we examine the flavor composition and angular dependence of QCD scattering 
channels (proportional to a 2 ). In the relevant region of Q, QCD jet production is dominated 
by quark-gluon scattering, qg — > qg, and, at lower Q values (not shown) by gluon-gluon 
scattering gg —> gg (cf. Fig. E](a)). Both types of contributions are peaked strongly at 
large \Ay\; see Fig. [2](b). The quark-quark scattering QCD contributions, qq — > qq and 
qq' — > qq 1 , have a flatter dependence on Q and \Ay\ than the gluon-scattering contributions. 
Similarly, the electroweak 0(a 2 EW + a s otEw) dijet production, also proceeding via quark- 
quark scattering, has a flatter dependence on \Ay\ than the gluon-scattering-dominated 
QCD background. 

The qq — > qq and qq' — > qq' QCD amplitudes can interfere with the electroweak am- 
plitudes having the same initial and final quark states. The largest contributions to the 
interference, together with the full interference term, are shown in Fig. [3] as a function of Q. 
The magnitude of the G(a s aEw) interference is similar to that of the 0(a 2 EW ) cross section. 
Its shape is determined by the interplay of the partial interference terms. 

We classify these partial terms according to their initial-state quark composition. That 
is, the cross section "uu, dd n corresponds to the sum of the interference cross sections with 
two initial-state u quarks or two d quarks, i.e., uu — > uu and dd — > dd\ "ud n corresponds 



3 Our results for an asymmetric cut —2 < yj < 1 on jet rapidities (exactly corresponding to the STAR 
acceptance) are qualitatively similar. 
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pp -> jjX, \[S = 500 GeV, QCD channels only 



pp -> jjX, \|S = 500 GeV 




Figure 2: 0{ot 2 s ) cross sections for pp — > jjX at = 500 GeV in various scattering channels 
(QCD contributions only), plotted versus (a) dijet invariant mass Q, (b) difference of jet rapidities 
Ay = \yi - yi\. 
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Figure 3: The unpolarized interference cross section (solid) and contributions of individual initial- 
state channels: uu and dd (short-dashed); ud (dotted); ud and du (dot-dashed); uu and dd (long- 
dashed) . 
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to the interference cross section for ud — > ud, ud — > du, du — > ud, and du — > du; and so 
on. Smaller interference terms involving strange or charm quarks are included in the total 
interference, but not shown separately. 

It is useful to distinguish between the interference terms that possess resonant properties 
and those that do not. The resonant interference occurs between an s-channel W or Z 
boson amplitude and a different conjugate amplitude with the same external states. The real 
part of the heavy boson propagator in the s-channel amplitude changes sign at the boson's 
mass pole, Q = My. Consequently a resonant interference term exhibits a pronounced 
enhancement slightly below the W or Z pole and a comparable enhancement of the opposite 
sign immediately above the pole. In our case, the "ud, du" cross section (dot-dashed curve) 
is the most prominent resonant term dominated by the interference with the s— channel W 
boson amplitude. It changes sign at Q = Mw- Despite its large magnitude, its contribution 
to the integrated event rate in the signal region 70 < Q < 90 GeV is small, due to the 
cancellation between the contributions of the opposite sign from below and above Q = M w . 
Similarly, a resonance driven by the s-channel Z boson amplitude occurs in the u uu, dd" 
cross section (long-dashed line) at Q = Mz- 

The largest interference contributions to the integrated rate in the signal region arise 
from non-resonant uu, ud, and dd processes (short-dashed and dotted lines). These cross 
sections do not contain a large resonant propagator, but rather are enhanced by a product 
of two large valence-dominated quark PDFs, u(x\, Q)u(x 2 , Q), etc. From this discussion, we 
conclude that the ud and du interference contributions cancel when integrated over the signal 
region. The largest surviving interference is due to the uu, ud, and dd terms. As a result 
the interference strongly affects dependence on u and d quark PDFs, but only marginally 
on the d antiquark PDF. 

C. Parton flavor composition of da/dy 

In Fig. S](a), we show partial contributions to the O (a EW ) cross sections da/dy at 70 < 
Q < 90 GeV involving only W + and W~ contributions in all channels. In this figure, 
we identify contributions proportional to u{x\), d{x\), u(x\), and d(xi). The pure W ± 
contribution is dominated by the resonant terms. Hence the curves in Fig.[4](a) closely follow 
the rapidity dependence suggested by Eqs. ([6]) and ([7]). Only one term [proportional to u(x) 
or d(x) at x — > 1] survives on the right-hand sides of these equations when y approaches its 
kinematical limits, y — > ±ha.(Q/y/s). 

The cross section supplied by W~ contributions constitutes about 1/3 of that from W + 
contributions. Consequently, the combined W ± cross section is dominated by u(x\) and 
d(xi) contributions at large positive and negative y, as in resonant W + production. 

Figure 01(b) shows the quark flavor decomposition for the electroweak+interference 
0{o? EW + a s a.Ew) cross section. In this case, the combined u(xi) contribution is roughly 
equal to the sum of the W + contribution in Fig. [4](a) and uu, ud interference contribution 
(symmetric with respect to y — 0). The combined d(x\) contribution is similarly made of 
the W~ and dd, du interference contributions. Both u(xi) and d(xi) contributions remain 
dominated by the resonant W terms, as the interference involving sea quarks is quite small. 
The essential conclusion to be drawn is that the combined electroweak+interference cross 
section largely preserves the PDF dependence of the resonant W + production, notably, the 
sensitivity to u{x\,Q) (d(xi,Q)) at the forward (backward) y values. 
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Figure 4: Quark flavor composition of the unpolarized cross section da/dy: (a) W{a 2 EW ) cross 
section only; (b) full electroweak+interference 0{a 2 EW + a s oiEw) cross section 



III. THE SINGLE-SPIN ASYMMETRY 

In Fig. [5] we present the single-spin asymmetry Ai(y), evaluated as a function of rapidity 
for two sets of DNS'2005 leading-order polarized PDFs and MRST'2002 NLO unpolarized 
PDFs. The statistical uncertainty 5A L in the measurement of the asymmetry is evaluated 
as 



SAi 



2 

beam 



(10) 



for the integrated luminosity C = 160 pb _1 and proton beam polarization Pbeam = 0.7. 

The asymmetry Al and its projected statistical uncertainty can be determined in two 
ways. In the first (most straightforward) method, the spin-averaged cross section is taken 
to be the full O (ct| + a s aEw + a Ew) cross section. The resulting rapidity (y) dependence 
of A L in the signal region 70 < Q < 90 GeV is shown in Fig. [5](a). The value of A L is 
generally small, and the statistical uncertainties (evaluated based on the total number of 
events for both beam polarizations in each bin) do not allow discrimination between the sets 
of polarized PDFs shown. 

In the second method, we set the unpolarized cross section equal to the O (a s a.Ew + &%w) 
cross section both in the denominator of Al and SAl. This procedure is roughly equivalent 
to the measurement of A L in which the background cross section (dominated by the 0(a 2 ) 
term) is measured precisely in the side bands (Q < 70 GeV and Q > 90 GeV), extrapolated 
into the signal region (70 < Q < 90 GeV), and subtracted from the measured unpolarized 
rate. Using this approximate side-band subtraction technique in Fig. [5](b), we see that the 
predicted magnitude of A L is increased substantially, and different parametrizations of the 
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Figure 5: Parity-violating asymmetry Al in the dijet mass interval 70 < Q < 90 GeV without 
and with side-band subtraction of the parity-conserving QCD background, plotted for two sets of 
leading-order DNS'2005 polarized PDFs. 



polarized PDFs can be discriminated, given the projected statistical uncertainties. 

The asymmetry Al includes both the W + and W~ contributions, as a result of summation 
over all jet charges. At large positive y, Al approximately follows —Au(x)/u(x) at large x, 
expected to satisfy Au(x)/u(x) — > 1 in the x — ► 1 limit. At large negative y, Al mostly 
reflects the behavior of sea quark PDFs at x — > 0.02, notably of Ad(x)/d(x). 

We now examine the region of negative y in greater detail in order to gauge the sensitivity 
of Al to variations in the polarized PDFs within the limits tolerated by current parametriza- 
tions. We plot the values of A L obtained after side-band subtraction for the DNS'05 NLO 
PDF sets with varied first moments of Au(x), Ad(x), Au(x), and Ad(x). The results are 
shown in Fig. El When y < —1, A L is clearly sensitive to the variation in the first moment of 
Ad and, to a smaller degree, of Am. It is nearly insensitive to the variations in Am and Ad. 
This result suggests that the measurement of Al will constrain Ad(x), given the projected 
statistical uncertainties. 

We remark that the pronounced sensitivity of Al to —Au(x)/u(x) and Ad(x)/d(x) arises 
because of the dominance of the W + contribution to parity- violating dijet production, nearly 
independently of the PDF set chosen. Equations (jH [Ml]) state that the resonant parity- 
violating part of Al behaves at large \y\ approximately as 

MeiI _l Mxi) r ( x , Xo ) 
lim x A L (y) - + ^ T ^ X2 \ (11) 



and 



■hx(M w /-Si) l+r(xi,x 2 ) 



l im A L (y)* ^», + ■*»' r( *" Jl) , (12) 



y->\n(M w /y/i) 1 + r(x 2 , Xi) 
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Figure 6: The asymmetry Al for backward production of jet pairs (y < 0), computed for DNS'2005 
NLO polarized PDFs [25|| after side-band subtraction. The black solid curve corresponds to the 
best-fit PDF set 1. The pairs of other curves correspond to the maximal and minimal values of 
Al obtained if the first moment Aq = f dxAq(x,Q = 3.16 GeV) deviates from its best-fit value 
by an amount corresponding to x 2 = 1.02x^ in , where Xmin ^ s the minimum of the log-likelihood 
function x 2 m the fit. We show Al for varied Ad (red long dashes); varied Au (brown short dashes); 
varied Au (blue dots); and varied Ad (green dot-dashes). At large negative rapidities, the most 
pronounced variations in Al are due to the variation of Ad. 



where 

r{xuX2)s d JfApl. (13) 
u{xi) d(x 2 ) 

The factor r is small (of order 1/3) in both limits, as a consequence of the smallness of the 
d(x)/u(x) ratio at x — > 1 (see, e.g., Figs. 3-6 in Ref. [H|)- Unless Au(xi)/u(xi) is much 
smaller than Ad{x\) / d{x\) in absolute magnitude, it dominates the numerator of A L (y) in 
Eq. (fTTl) . Similarly, Ad(xx) / d{xi) dominates the numerator of Eq. (Tl2l) . unless it is much 
smaller than Au(xi)/u(xi) in absolute magnitude. 

In the lepton decay mode, the positron from a W + decay tends to scatter into the central 
rapidity region, hence smearing the x dependence of the underlying d parton distribution 
14 30l |. In jet production, the smearing does not occur as a result of the direct measurement 



of the rapidity y of the jet pair. 

Observation of Al in the jet pair mode, after the background subtraction, appears to 
guarantee a large asymmetry associated with Au(x) at x — > 1 for forward jet pair rapidities 
and to test Ad(x) in the region of negative jet pair rapidities. 
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IV. CONCLUSIONS AND DISCUSSION 



Data on W production obtained from longitudinally polarized proton-proton scattering 
at RHIC can be used to extract spin-dependent quark and antiquark parton distributions 
Aqi(x, M w ) and Aq^x, M w ), providing information complementary to that obtained from 
polarized deep-inelastic lepton scattering. The longitudinal spin asymmetry A L (y) of the 
rapidity (y) dependence of W production is particularly sensitive to the spin-dependent 
quark PDFs. In this paper, we examine the prospects for definitive measurements of Ai(y) 
when the W boson is detected at RHIC in its hadronic decay mode, W — > 2 jets, at the 
collision energy y/s = 500 GeV. 

To take advantage of the increase in the event rate and, in principle, of the direct mea- 
surement of y offered by the hadronic decay mode, one must first address the challenges of 
the large parity-conserving background from QCD and electromagnetic production of jets 
in the vicinity of the W signal. In Sec. [Til of this paper, we report our tree-level computa- 
tion of the amplitudes for all electroweak and strong interactions processes that lead to a 
pair of jets in proton-proton scattering. We treat carefully all effects of interference among 
the subprocesses. As shown in Fig. CD, the backgrounds in unpolarized scattering appear 
admittedly daunting. Most of this background is due to subprocesses with emission of one 
or two final-state gluons, as shown in Fig. El Nevertheless, even at the unpolarized level, an 
experimental side-band subtraction procedure could be used to reach acceptable sensitivity 
to events containing two final-state quarks from electroweak scattering or QCD-electroweak 
interference. With this procedure in place, the unpolarized event rate is dominated by pro- 
cesses involving the u quark PDF in the forward region of the dijet y, and by processes 
involving the d antiquark PDF in the backward region of y, even without separation of 
dijets from W + and W~ production [cf. Fig. Hl(b)]. 

Turning to the longitudinally polarized case, we remark first that the largest QCD and 
electromagnetic backgrounds cancel in the parity- violating numerator cr(p~*p) — cr(p^p) of 
A L (y). This means that the desired sensitivity of A L (y) to Au(x, M w ) and Ad(x, M w ) is 
mostly preserved. The magnitude of A L (y) can be enhanced by applying the above sub- 
traction procedure to the spin-averaged denominator a(p~*p) + cr(p*~p) of Ai(y). The strong 
and electromagnetic backgrounds may be measured in the dijet invariant mass distribu- 
tion in side-band regions next to the W signal location. An appropriate extrapolation and 
subtraction of the backgrounds may then be performed when defining the denominator of 
A L (y). 

Our quantitative predictions for the asymmetry are presented in Figs. [5]and[6l In Fig. [5], 
we demonstrate that, once the side-band subtraction technique is employed, the predicted 
magnitude of Al is sufficiently large, when compared to the expected uncertainties. The 
data on A L (y) would thus be a very valuable component of a global analysis leading to 
definitive longitudinally polarized PDFs. 

At large positive y, Al approximately tracks —Au(x)/u(x) at large x, commonly assumed 
to satisfy Au(x)/u(x) — > 1 at x — > 1. We find that Al can drop as low as —0.5 at y — 1.2—1.5 
(Fig. El(b)). The statistical error in this range ofy is still small enough to allow a meaningful 
measurement. 

At large negative y, A L reflects the behavior of the sea quark density Ad(x)/d(x). One 
way to estimate the power of discrimination is to compute the changes in the predicted 
Al(v) that result from varying the polarized parton densities within the bands of tolerance 
of current parametrizations. As shown in Fig. EH the most pronounced predicted variations 
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in the backward rapidity region result from changes in Ad. 

In an effort to increase the magnitude of Ai(y), several techniques may prove to be 
helpful. Selection of events at central Ay = \yi — 2/ 2 1 (or central cos#* in the dijet rest 
frame) helps to suppress the gluon-dominated QCD background, as follows from Figs. [T](c) 
and [2(b). 



One may consider selecting jets containing a final-state charmed particle [24|. In this 
case, the QCD background is reduced much more than the W signal rate by excluding many 
sources of light-flavored jets. Using realistic efficiencies for charm detection, we cannot 
conclude that the charm tag clearly improves the expected significance of the measurements 
of Ai(y). An alternative approach could involve a tag on a leading charged particle in 
the jet (e.g., a charged pion carrying at least a third of the total jet energy) to emphasize 
contributions of the quark-initiated jets over gluon-initiated jets. As with the charm tags, the 
reduction in the background over signal in this case must be balanced against the substantial 
suppression in the total event rate. We therefore defer our conclusion about the feasibility 
of charm and leading-particle tags until a more detailed study. 

The analysis reported here uses hard-scattering cross sections computed at the lowest 
order in perturbation theory. Next-to-leading order (NLO) contributions in QCD increase 
the predicted rates for hadronic jet production 3ll, [32j and are essential for reducing large 
scale dependence of the lowest-order cross section. This increase means that the backgrounds 
in the denominator of Ajiy) can be larger than our estimates. On the other hand, as 
shown by Moretti et al. [33fl, comparably large enhancements are predicted in the single spin 
polarized cross sections at RHIC, meaning naively that the predicted magnitude of Al(?/) 
could remain largely unaffected. A quantitative statement about the ultimate effects of NLO 
terms on Ai(y) would require a consistent treatment of NLO contributions to all contributing 
processes, in both the unpolarized and polarized well as consistent application of 

kinematic selections that we find advantageous. We leave this task for future work, but we 
note here that comparably large enhancements of both the polarized and unpolarized rates 
would lead to a decrease in the projected statistical uncertainty SA L , Eq. (fTUl) . since 5A L is 
proportional to 1/yN, where N is the number of events observed. 

Our essential conclusion is that an experimental study of high-pr jet pair production at 
RHIC appears to offer a promising way to make definitive measurements of the longitudinal 
spin asymmetry Ai(y). 
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